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© The characteristics of a biological fluid sample 
having an analyte are determined from a model 
constructed from plural known biological fluid sam- 
ples. The model is a function of the concentration of 
materials in the known fluid samples as a function of 
absorption of wideband infrared energy. The 
5! wideband infrared energy is coupled to the analyte 
^containing sample so there is differential absorption 
Wof the infrared energy as a function of the 
JO wavelength of the wideband infrared energy incident 
on the analyte containing sample. The differential 
absorption causes intensity variations of the infrared 
® energy incident on the analyte containing sample as 
a function of sample wavelength of the energy, and 
©concentration of the unknown analyte is determined 
^from the thus-derived intensity variations of the in- 
UJfrared energy as a function of wavelength from the 
model absorption versus wavelength function. 
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METHOD OF AND APPARATUS FOR DETERMINING THE SIMILARITY OF A BIOLOGICAL ANALYTE FROM 
A MODEL CONSTRUCTED FROM KNOWN BIOLOGICAL FLUIDS 



Field of Invention 

The present invention relates generally to de- 
termining the nature, i.e., the similarity or con- 
centration, of a biological analyte in comparison 
with a model constructed from plural known bio- 
logical fluids, and, more particularly, to such a 
method and apparatus wherein a sample of biologi- 
cal fluid containing the analyte is irradiated with 
infrared energy having at least several wavelengths 
to cause differential absorption by the analyte as a 
function of the wavelengths and properties of the 
analyte. 



Background Art 

For various care and treatment of mammal 
patients, it is necessary to determine concentra- 
tions of certain species in biological fluids. For 
instance, diabetics must be apprised of their blood 
glucose concentrations to enable insulin dosage to 
be adjusted. To determine blood glucose concen- 
trations, blood is presently drawn several times per 
day by the diabetic, usually via a finger prick. If the 
blood glucose concentrations in such individuals 
are not properly maintained, the individuals be- 
come susceptible to numerous physiological prob- 
lems, such as blindness, circulatory disorders, 
coronary artery disease, and renal failure. For 
these reasons, a substantial improvement in the 
quality of life of persons suffering from various 
maladies, such as diabetes meliitus, could be at- 
tained if the concentrations of species in body 
fluids are non-invasively and/or continuously deter- 
mined. For example, for diabetic patients having 
external or implantable insulin pumps, a feedback 
loop for these pumps could be controlled by con- 
tinuously monitoring glucose concentrations, to en- 
able an artificial pancreas to be developed. 

Exemplary systems have been previously pro- 
posed to monitor glucose in blood, as is necessary, 
for example, to control diabetic patients. This prior 
art is represented, for example, by Kaiser, U.S. 
Patent 4,169.676, Muiler, U.S. Patent 4,427,889, 
and Dahne et al, European Patent Application 
84810212.5, and Bauer et al, Analytica Chimlca 
Acta 197 (1987) pp. 295-301. 

In Kaiser, glucose in blood is determined by 
irradiating a sample of the blood with a carbon 
dioxide laser source emitting a coherent beam, at a 
single frequency, in the mid-infrared region. An 
infrared beam derived from the laser source is 
coupled to the sample by way of an attenuated 



total reflectance crystal for the purpose of contact- 
ing the blood sample. The apparatus uses double 
beam instrumentation to examine the difference in 
absorption at the single frequency in the presence 

5 and absence of a sample. The reliability of the 
Kaiser device is materially impaired in certain situ- 
ations because of the reliance on a single fre- 
quency beam for reasons explained below. Also, 
we have found from calculations based on available 

w information that Kaiser's statement anent optical 
energy penetrating the skin to the depth of the 
blood capillaries is unlikely due to water absorption 
of the mid-infrared beam. 

Muller discloses a system for quantifying glu- 

75 cose in blood by irradiating a sample of the blood 
with energy in a single beam from a laser operating 
at two frequencies in the mid-infrared region. The 
infrared radiation is either transmitted directly to 
the sample or by way of an attenuated total reflec- 

20 tance crystal for in vitro sampling. One frequency 
that irradiates the sample is in the 10.53 - 10.65 
micrometer range, while the other irradiating fre- 
quency is in the 9.13-9.17 micrometer range. The 
radiation at the first frequency establishes a 

25 baseline absorption by the sample, while glucose 
absorption by the sample is determined from the 
intensity reduction caused by the sample at the 
second wavelength. The absorption ratio by the 
sample at the first and second frequencies quanti- 
se? fies the glucose of the sample. There is no glucose 
absorption at the first wavelength. 

Dahne et al employs near-infrared spectro- 
scopy for non-invasively transmitting optical energy 
in the near-infrared spectrum through a finger or 

35 earlobe of a subject. Also discussed is the use of 
near-infrared energy diffusely reflected from deep 
within the tissue. Responses are derived at two 
different wavelengths to quantify glucose in the 
subject. One of the wavelengths is used to deter- 

40 mine background absorption, while the other 
wavelength is used to determine glucose absorp- 
tion. The ratio of the derived intensity at the two 
different wavelengths determines the quantity of 
glucose in the analyte biological fluid sample. 

45 Bauer et al discloses monitoring glucose 
through the use of Fourier-transform infrared spec- 
trometry wherein several absorbance versus 
wavelength curves are illustrated. A glucose con- 
centration versus absorbance calibration curve, dis- 

50 cussed in the last paragraph on p. 298, is con- 
structed from several samples having known con- 
centrations, in response to the intensity of the 
infrared energy absorbed by the samples at one 
wavelength, indicated as preferably 1035 cm" 1 . 
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All of the foregoing prior art techniques thus 
use only a single frequency analysis or ratio of two 
frequencies to determine a single proportionality 
constant describing a relationship between absorp- 
tion of the infrared energy by the sample and 
concentration of a constituent of the biological fluid 
sample being analyzed, usually glucose. Hence, 
the prior art analysis is univariate since absorption 
by the constituent of interest at a single wavelength 
is determined. 

However, univariate analysis has a tendency to 
be inaccurate in situations wherein there are con- 
centration variations of any substance which ab- 
sorbs at the analysis frequency. Biological systems 
are subject to numerous physiological perturbations 
over time and from person to person. The perturba- 
tions cause inaccuracies in univariate analysis, 
thereby decreasing the accuracy and precision of 
such analysis. The physiological perturbations in- 
volving any substance which absorbs at the analy- 
sis frequencies do not permit an operator of a 
system utilizing univariate analysis to recognize the 
resulting inaccuracy. In addition, nonlinearities may 
arise from spectroscopic instrumentation, refractive 
index dispersion, or interactions between molecules 
of the sample which cannot generally be modelled 
by univariate techniques. In addition, unknown bio- 
logical materials in the sample have a tendency to 
interfere with the analysis process, particularly 
when these materials are present in varying 
amounts. Also the univariate techniques are usually 
not capable of identifying outlier samples, i.e.. sam- 
ples with data or constituents or spectra among the 
calibration or unknown data which differ from the 
remainder of the calibration set 

The described prior art systems utilizing mid- 
infrared energy are not feasible for non-invasive in 
vivo determinations of glucose concentrations b£ 
cause of penetration depth limitations. 

The most frequently employed prior art tech- 
niques for determining the concentration of molec- 
ular substances in biological fluids have used en- 
zymatic, chemical and/or immunological methods. 
However, all of these techniques require invasive 
methods to draw a blood sample from a subject; 
typically, blood must be drawn several times a day 
by a finger prick, such as presently employed by a 
diabetic. For example, in the determination of glu- 
cose by diabetics, such invasive techniques must 
be performed using present technology. It would 
be highly desirable to provide a less-invasive, con- 
tinuous or semi-continuous system for automati- 
cally analyzing glucose concentrations in the con- 
trol of diabetes mellitus. 

It is, accordingly, an object of the present in- 
vention to provide a new and improved method of 
and apparatus for determining characteristics of a 
biological analyte sample. 



Another object of the present invention is to 
provide a new and improved apparatus for and 
method of using infrared energy for analyzing bio- 
logical fluids wherein the apparatus and method 

5 are particularly suitable for analyzing samples hav- 
ing concentrations of substances which variably or 
differentially absorb the infrared energy. 

Another object of the invention is to provide a 
new and improved method of and apparatus for 

10 utilizing infrared energy to determine a characteris- 
tic, e.g., concentration, of a biological analyte by 
comparison of the absorption characteristics of said 
sample with a mathematical model constructed 
from several pre-stored biological fluids having 

75 known absorption versus wavelength characteristics 
at known analyte concentrations. 

A further object of the invention is to provide a 
new and improved apparatus for and method of 
analyzing biological fluids with infrared energy 

20 wherein interference with the infrared energy due 
to numerous physiological perturbations over time 
and between people does not have a particularly 
adverse effect on the results. 

An additional object of the invention is to pro- 

25 vide a new and improved apparatus for and meth- 
od of using infrared energy to analyze biological 
fluids, wherein non-linearities due to various 
causes, for example, spectroscopic instrumenta- 
tion, refractive index dispersion, and/or inter-molec- 

30 ular interactions, do not have an adverse effect on 
the analysis results. 

An additional object of the present invention is 
to provide a new and improved apparatus for and 
method of using infrared energy to determine the 

35 nature of a biological sample wherein the presence 
of unknown biological materials in the sample does 
not interfere with the analysis of the sample, as 
long as these unknown biological materials are 
present in a calibration set which is used to derive 

40 a mathematical model which represents the re- 
sponse of known fluids to the infrared energy. 

A further object of the invention is to provide a 
new and improved apparatus for and method of 
using infrared energy to determine characteristics 

45 of biological fluids wherein outlier samples subsist- 
ing in a calibration set used to derive a mathemat- 
ical model are identified and either eliminated or 
accommodated so as not to have an adverse effect 
on the determination. 

so Another object of the invention is to provide a 
method of and apparatus for identifying the pres- 
ence of outliers. The quality of the calibration re- 
sults and the reliability of the unknown sample 
analyses can be critically dependent on the detec- 

55 tion of outlier samples. In the calibration set, an 
outlier is a sample that does not exhibit the char- 
acteristic relationship between composition and the 
absorbance spectrum of the other calibration sam- 
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pies. During prediction, an outlier is a sample that 
is not representative of samples in the calibration 
set. Outliers in the calibration samples can impair 
the precision and accuracy of the calibration and 
limit the quality of the analyses of all future sam- 
ples. The results of the analyses of outlier unknown 
samples by multivariate calibration cannot be con- 
sidered reliable, and samples containing outliers 
should be analyzed by other methods. Thus, effi- 
cient detection of outlier samples is crucial for the 
successful application and wide acceptance of mul- 
tivariate spectral analyses. For example, outliers 
occur as a result of changes in instrumental re- 
sponse, incorrect analyte determination by the ref- 
erence method, unique type off sample, unexpect- 
ed components, unusual baseline, incorrectly label- 
ed or documented sample, etc. 

Still an additional object of the invention is to 
provide a new and improved biological fluid analy- 
sis apparatus and method which is particularly 
adaptable, in certain embodiments, to non-invasive 
determinations. 



The invention 

In accordance with one aspect of the present 
invention, the concentration of a biological fluid 
containing an analyte is determined from a model 
constructed from plural known biological fluid sam- 
ples. The model is a function of the concentration 
of materials in the known samples as a function of 
absorption at at least several wavelengths of in- 
frared energy. The infrared energy is coupled to 
the analyte containing sample so there is differing 
absorption of the infrared energy as a function of 
the several wavelengths and characteristics of the 
analyte containing sample. The differing absorption 
causes intensity variations of the infrared energy 
passing through the analyte containing sample as a 
function of the several wavelengths. The thus-de- 
rived intensity variations of the infrared energy as a 
function of the several wavelengths are compared 
with the calibration model relating concentration to 
plural absorption versus wavelength patterns de- 
rived from the plural known biological fluid samples 
having various concentrations. The comparison en- 
ables the determination of the analyte concentra- 
tion from the measured intensity variations for the 
biological fluid containing the analyte. In the pre- 
ferred embodiment, the comparison is made in a 
computer by the partial least squares method, al- 
though other multivariate techniques can be em- 
ployed. 

In the computer implementation, the intensity 
variations as a function of wavelength are con- 
verted into plural first electric signals, such that 
different ones of the first electric signals are as- 



signed to different ones of the wavelengths. The 
magnitude of each of the different first signals is 
determined by the intensity of the transmitted en- 
ergy at the wavelength assigned to that particular 

5 first signal. The transmitted energy in the presence 
of the analyte containing sample is statistically 
compared with the transmitted energy in- the ab- 
sence of the sample to determine the 'absorption 
by the biological analyte containing fluid. 

70 A multivariate statistical method, preferably us- 
ing the partial least squares technique in a manner 
known in the statistical art, enables a model to be 
constructed of the infrared absorption versus 
wavelength characteristics and analyte concentra- 

75 tions. Following determination of the calibration 
model, the infrared absorption versus wavelength 
of the unknown fluid enables estimation of the 
analyte concentration. For example, if blood is be- 
ing monitored, the glucose concentration of the 

20 blood can be determined by a statistical compari- 
son between the model and the unknown sample 
absorption of the infrared energy at several 
wavelengths. If the absorption versus wavelength 
characteristics, i.e., spectrum, of the unknown fluid 

25 containing the analyte differ sufficiently from the 
spectrum generated by the model, the statistical 
analysis enables a determination to be made that 
the absorption of the unknown fluid at the different 
wavelengths does not closely enough match the 

30 model to provide any meaningful data. This is 
particularly important for certain applications, e.g., 
relating to insulin control for a subject. A controller 
for insulin injection is not activated in response to 
the absorption versus wavelength characteristics 

35 derived from the analyte if the analyte data do not 
adequately fit the model. 

It has been found that the presence of alcohol 
in the blood of a diabetic overlaps with the analyte 
absorption versus wavelength characteristics to 

40 such an extent that accurate determination of glu- 
cose concentration is impossible with univariate 
methods, as in the prior art. By utilizing the mul- 
tivariate technique of the present invention, it is 
possible to derive an indication that the absorption 

45 data for the unknown sample at several frequencies 
do not provide a basis upon which to derive mean- 
ingful glucose concentrations for determination of 
control signals to an insulin pump. Thereby. If the 
measured unknown spectrum differs to such an 

so extent from the model that there is, in fact, no 
adequate fit between them, the insulin pump con- 
troller is not modified and the operator is made 
aware of a possible error in the concentration de- 
termination. If the partial least squares method is 

55 used to determine the fit between the model and 
the unknown sample absorption versus wavelength 
response, a statistical measure representing the 
quality of the fit between the model and the un- 
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known spectrum in excess of a predetermined val- 
ue provides an indication that analysis may be 
unreliable. 

In accordance with a further aspect of the 
invention, outlier samples can be detected and/or 
eliminated from a calibration set of biological fluids. 
The outlier samples are detected and eliminated 
from the calibration set by utilizing the multivariate 
technique for each calibration sample forming the 
model. The absorption versus wavelength re- 
sponses for at least several wavelengths of each of 
the fluids forming the calibration set are derived. 
The absorption versus wavelength characteristic of 
each sample is compared to the model generated 
from all of the other calibration samples. A pre- 
determined level of acceptable variation for each 
sample spectrum from the model is established. A 
statistical test (typically the calculation of F-statis- 
tics or the Mahalanobis distance) is employed to 
determine which sample or samples of the calibra- 
tion set exceed the predetermined level of accept- 
able variation. Hence, for example, if a particular 
fluid contemplated for use in establishing the cali- 
bration set differs from the remaining samples in 
the set by a predetermined amount, that particular 
fluid is not used as a member of the calibration set. 

The respective intensity variations can be con- 
verted into electric signals using several different 
frequency dispersive techniques including Fourier 
transform techniques. One possible realization uti- 
lizes an array of infrared-to-electrical transducers, 
one of which is provided for each of the 
wavelengths so that different ones of the trans- 
ducers derive different ones of the electric signals. 
In one of the embodiments, the array includes 
multiple individual filters each having a passband 
for one of the wavelengths. The filters are posi- 
tioned relative to the transducers so that the in- 
frared energy passed through one of the filters is 
incident on a corresponding one of the transducers. 
In a second embodiment, the array includes a 
sheet of infrared gradient wavelength responsive 
material having areas positioned relative to the 
transducers so that the infrared energy passed 
through different areas of the sheet is incident on 
corresponding different transducers. In a third em- 
bodiment, each of the transducers has a different 
chemical composition or is doped to be responsive 
to a different one of the several wavelengths. 

In one embodiment, particularly adapted for 
monitoring in vitro biological fluids, the infrared 
energy is coupled to the fluid via an attenuated 
total reflectance crystal and the infrared source is 
in the mid-infrared spectrum. It is also possible to 
use external mid-infrared sampling for biological 
fluids which are circulating outside of the body, 
such as blood components during dialysis. For in 
situ analysis, such a crystal is implanted in the 



body and has a bio-compatible coating saturated or 
in contact with body fluids. 

In a second embodiment, a catheter including 
a fiber optic element is inserted into an artery of a 

5 subject. The end of the fiber optic element in the 
artery has a reflective end coating. Infrared energy, 
either in the near-infrared or in the mid-infrared 
spectrum, illuminates the opposite end of the fiber 
optic element. The portion of the fiber optic ele- 

70 ment in the artery functions in the same manner as 
an attenuated total reflectance crystal. In systems 
utilizing the attenuated total reflectance element 
and the fiber optic element, blood or other body 
fluids irradiated by the optical energy of the in- 

15 frared source differentially absorb the different fre- 
quencies to enable the multivariate analysts to be 
performed. 

In a third embodiment, an optical fiber element 
transmitting in either the mid-infrared or near-in- 

20 frared spectrum is inserted into the body for analy- 
sis of interstitial or subdermal fluid. The portion of 
the fiber in contact with the interstitial fluid func- 
tions in the same manner as an attenuated total 
reflectance element. In the specifically disclosed 

25 embodiments, such techniques are used for fluid 
monitoring in an earlobe through the use of a fiber 
optic loop including a fiber optic element that 
pierces the earlobe and transmits the non-absorbed 
energy to a frequency dispersing array. 

30 In other embodiments of the invention, optical 
energy in the near-infrared spectrum is directly 
transmitted through a body portion and an analysis 
of wavelength versus absorption over several 
wavelengths of the spectrum is performed. In the 

35 specifically disclosed embodiments, a digit, such 
as the index finger of a subject, or an artery, is 
irradiated with several wavelengths in the near- 
infrared spectrum. 

In a further embodiment, a biological sample 

40 fluid in an internal organ, such as the brain or liver, 
is irradiated by several wavelengths in the near- 
infrared spectrum. A source of the near-infrared 
energy is placed in close proximity to a body part, 
such as the head or abdomen, containing biological 

45 fluid. The sample fluid in the organ diffusely re- 
flects the near-infrared wavelengths and intensity 
versus wavelength responses at several 
wavelengths are measured and used in with the 
derived mathematical model to derive the analyte 

50 concentration. 

Hence, the sample being analyzed includes a 
component that absorbs the infrared energy at a 
plurality of the wavelengths to which the source- 
detector combination is sensitive so that the thus- 

55 derived intensity variations are determined by ab- 
sorption by the component at the plural 
wavelengths. The plural pre-stored patterns repre- 
sent absorption by the known fluids at these plural 
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wavelengths. The concentration determination step 
involves comparing the intensity of the thus-derived 
intensity variations at the plural wavelengths with 
mathematical model derived from the absorption 
intensities of the pre-stored patterns at these plural 
wavelengths. The component, in a preferred em- 
bodiment, consists of glucose. However, the device 
could be used for quantizing alcohol, ketones, 
fatty acids, cholesterol, lipoproteins, triglycerides, 
white blood cells, albumin, blood urea nitrogen, 
creatinine, concurrent medications such as drugs, 
inorganic molecules such as phosphates, and any 
other infrared absorptive component present in a 
biological analyte fluid. 

Other source/detector configurations are fea- 
sible, such as a single or multiple element detector 
in conjunction with a dispersive Hadamard trans- 
former or a Fourier transform infrared spectrom- 
eter. In addition, a standard prism or grating as- 
sembly could be employed. In the disclosed em- 
bodiment, there is a frequency dispersing arrange- 
ment to provide a unique frequency or wavelength 
intensity indication at each array element or chan- 
nel. The disclosed structures for frequency disper- 
sal, however, are of particular value for the present 
application because they are relatively small 
dispersion/detector assemblies. 

The mid-infrared spectrum is preferred in ap- 
plications utilizing attenuated total reflectance cry- 
stals or fiber optic elements because the mid- 
infrared frequencies have a severely limited pene- 
tration depth into biological fluids. The penetration 
depth of the mid-infrared spectrum into biological 
fluids occurs primarily as a result of water absorp- 
tion so that transmission methods are not usually 
feasible in the mid-infrared spectrum. Hence, for 
embodiments of the invention using transmissivity, 
rather than reflectance, techniques, the infrared 
source is preferably in the near-infrared spectrum. 

Typically, absorptivities of biological fluids in 
the mid-infrared spectrum are about an order of 
magnitude more intense than in the near-infrared 
spectrum. This is because the spectral features in 
the near-infrared spectrum are overtones and com- 
bination bands of spectral features in the mid- 
infrared region. While attenuated total reflectance 
methods could be used with a near-infrared source, 
the decreased absorptivity of molecular species in 
the near-infrared spectrum decreases the precision 
of the analysis, without any benefits relative to 
using energy in the mid-infrared region. Hence, for 
embodiments of the invention using attenuated to- 
tal reflectance, the infrared source is preferably in 
the mid-infrared spectrum. 

Certain embodiments of the invention are par- 
ticularly adapted to be used for in situ investiga- 
tions, wherein a device is implanted in the body of 
a subject. The embodiment utilizing the mid-in- 



frared spectrum and the bio-compatible attenuating 
total reflectance crystal, the embodiment which 
monitors blood flow through an artery and the 
embodiment which monitors interstitial or subder- 

s mal fluids are designed to be implanted in the 
human body cavity. Such devices are particularly 
adapted to provide continuous control signals to an 
insulin pump of a diabetic. The insulin pump re- 
sponsive to a detector of the present invention 

70 replaces the physiological activities of the pancreas 
by simultaneously sensing glucose levels and pro- 
viding proper insulin delivery. 

The above and still further objects, features 
and advantages of the present invention will be- 

15 come apparent upon consideration of the following 
detailed description of a specific embodiment 
thereof, especially when taken in conjunction with 
the accompanying drawings. 

20 

Brief Description of Drawing 

Fig. 1 is a schematic diagram of a first 
embodiment of the invention, particularly adapted 

25 for monitoring in vitro biological fluids using mid- 
infrared spectroscopy and an attenuated total re- 
flectance crystal; 

Fig. 2 is a schematic diagram of another 
embodiment of the invention employing an im- 

30 planted sensor for monitoring biological fluids using 
mid-infrared spectroscopy and an attenuated total 
reflectance crystal; 

Rg. 3 is a schematic diagram of a further 
embodiment of the invention wherein blood in an 

35 artery is sampled via a fiber optic element for a 
mid-infrared or a near-infrared source; 

Rg. 4 is a schematic diagram of a further 
embodiment of the invention using a fiber optic 
loop for monitoring biological fluids in an earlobe of 

40 a subject, wherein the earlobe is irradiated via the 
fiber optic loop by a source emitting several 
wavelengths of optical energy in the near-infrared 
or mid-infrared spectrum; 

Fig. 5 is a schematic, perspective view of a 

45 further embodiment of the invention particularly 
adapted for monitoring blood components in a hu- 
man index finger using near-infrared spectroscopy 
and non-invasive transmission sampling; 

Fig. 6 is a schematic, perspective view of a 

so further embodiment of the invention for monitoring 
biological fluids in the brain non-invasively by using 
near-infrared spectroscopy and diffuse reflectance 
sampling; 

Fig. 7 is a schematic, perspective view of 
55 another embodiment of the invention wherein a 
body implanted housing containing a near-infrared 
source irradiates blood in an artery of a subject; 

Rg. 8 is a schematic, perspective view of a 
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first embodiment of the invention employing mul- 
tiple individual filter elements for providing fre- 
quency dispersion; 

Fig. 9 is a schematic, perspective view of a 
second embodiment of a frequency dispersion de- 
vice wherein a single gradient response filter sheet 
controls the wavelength of infrared energy incident 
on elements of a diode array detector; and 

Fig. 10 is a perspective view of a further 
embodiment of the invention wherein each element 
of a diode detector array is constructed so that it is 
responsive to optical energy of a different infrared 
wavelength. 

Description of the Preferred Embodiments 

Reference is now made to Fig. 1 of the draw- 
ing wherein vessel 1 1 contains human blood sam- 
ple 12 to be analyzed for glucose content. Blood 
analyte sample 12 is obtained from a typical in 
vitro source, such as from a patient undergoing 
surgery in a hospital operating room or a patient in 
a hospital intensive care ward. Blood analyte flows 
into and out of vessel 11 via conduits 8 and 9. 
respectively. Extending between and through end 
wails of vessel 11 is attenuating total reflectance 
crystal 13, having an index of refraction greater 
than the index of refraction of blood sample 12. 
Crystal 13 includes conical, opposite end portions 
14 and 15 which extend beyond end walls of 
vessel 11. O-ring seals 16 secure crystal 13 in 
place in apertures of the end walls of vessel 11, to 
prevent leakage of analyte sample 12 from the 
vessel. 

Optical energy from wideband mid-infrared 
(between 50 and 2.5 micrometers) black body 
source 17 is coupled to tapered wall portion 14 via 
cylindrical lens assembly 18. Lens assembly 18 
includes central convex reflector surface 19, having 
an axis coincident with the longitudinal axis of 
crystal 13. Convex surface 19 faces source 17, that 
is aligned with the axis of crystal 13. Lens assem- 
bly 18 comprises a circular mirror including con- 
cave reflecting surfaces 21 and 22, respectively 
positioned above and below the longitudinal axis of 
crystaJ 13. Concave reflecting surfaces 21 and 22 
are disposed on opposite sides, above and below, 
reflecting surface 19. Reflecting surfaces 19, 21 
and 22 are such that optical energy from source 17 
is initially incident on reflecting surface 19, thence 
is reflected to concave reflecting surfaces 21 and 
22 and thence is directed to conical end portion 14 
of crystal 13. The optical energy from source 17 
coupled to crystal 13 is reflected in the crystal as 
shown by ray paths 23 and 24. The structure 
illustrated in Rgure 1 including vessel 11, crystal 
13 and lens assembly 18 is commercially available 
from Spectra-Tech, Inc., under the CIRCLE CELL 



trademark. 

Some of the optical energy coupled from 
source 17 to crystal is transmitted from the crystal 
to sample 12, as shown by wavy lines 25 just 

5 beyond the intersection of the crystal and sample. 
Different wavelengths, i.e.. frequencies, of source 
17 are differentially absorbed by sample 12 as a 
function of concentrations of organic and inorganic 
molecules in the sample. In other words, certain 

io wavelengths emitted by source 17 are absorbed by 
sample 12 to a greater extent than other 
wavelengths from the source. The absorption oc- 
curs in sample 12 just beyond the interface be- 
tween crystal 13 and the sample. Since the degree 

J5 of absorption determines the amount of reflectance, 
those wavelengths which are most highly absorbed 
have the lowest intensity as they propagate from 
tapered, conical end portion 15 of crystal 13. 

The intensity versus wavelength characteristics 

20 of the infrared energy for at least several of the 
wavelengths of source 17 propagated from end 
portion 14 of crystal 13 are detected to enable the 
concentration of sample 12 to be determined. To 
this end, the infrared energy propagating from 

25 conical end face portion 15 of crystal 13 is incident 
on concave reflecting surfaces 31 and 32 of lens 
18. From reflecting surfaces 31 and 32, the non- 
absorbed infrared energy from source 17 is in- 
cident on the convex reflecting surface 33. From 

30 reflecting face 33, the optical energy is coupled to 
a frequency dispersion device 34 including diode 
array detector 35. 

Frequency dispersion device 34 and 
photodiode array detector 35 are arranged so that 

35 the array includes multiple output leads, one of 
which is assigned to a particular wavelength or 
narrow range of wavelengths of source 17. The 
amplitude of the voltage developed on each of the 
leads is commensurate with the intensity of the 

40 infrared energy incident on each particular detector 
in array 35 for the wavelength of the source asso- 
ciated with that detector. Typically, the photodiodes 
of array detector 35 are passive, rather than 
photovoltaic although photovoltaic devices may be 

45 employed. Thereby, the diodes of array detector 
35 must be supplied with DC power supply vol- 
tages, as derived from power supply 36 and coup- 
led to the diodes of array detector 35 via cable 37. 
The impedance of the diode elements of array 

50 detector is changed as a function of the intensity of 
the optical energy incident thereon in the passband 
of source 17 associated with each particular 
photodiode element. The impedance changes con- 
trol the amplitude of signals supplied by array 

55 detector 35 via bus 38 to a random access mem- 
ory of computer 39. 

Computer 39 includes a memory having stored 
therein a multivariate statistical model determined 
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from the concentration versus absorbance versus 
wavelength characteristics of several known blood 
samples. Such a model is constructed using tech- 
niques known by statisticians. 

Computer 39 determines the biological char- 
acteristics of analyte sample 12 by comparing the 
model with the amplitude versus wavelength re- 
sponse from array 35. Preferably, the comparison 
is made by a partial least squares technique, as 
disclosed, for example, by Lindberg et al, Analyt- 
ical Chemistry , Vol. LV, p. 643 (1983) in an article 
entitled, "Partial Least-Squares Method for Spec- 
trofluorometric Analysis of Mixtures of Humic Acid 
and Ligninsulfonate;" Martens et al, Applied Spec- 
troscopy , Vol. XL, p. 303 (1986) entitled "Near- 
infrared Reflectance Determination of Sensory 
Quality of Peas;" Haaland et al, Analytical Chem- 
istry , Vol. LX, p. 1193 (1988) in an article entitled 
"Partial Least Squares for Spectra! Analyses. 1. 
Relation to Other Quantitative Calibration Methods 
and the Extraction of Qualitative Information;" Haa- 
land et al, Analytical Chemistry , Vol. LX, p. 1203 
(1988) in an article entitled "Partial Least-Squares 
Methods for Spectral Analyses. 2. Application to 
Simulated and Glass Spectral Data;" Haaland, Ana- 
lytical Chemistry , Vol. LX, p. 1208 (1988), in an 
article entitled "Quantitative Infrared Analysis of 
Brophosphosilicate Films Using Multivariate Statis- 
tical Methods;" and Haaland et al, Proceedings of 
Pittsburgh Conference , Vol. XL, p. 874 (1989) in an 
article entitled "Outlier Detection During Multi- 
variate Quantitative Analysis of Spectroscopic 
Data." While the partial least squares method has 
been found to be precise, other techniques, such 
as a principal component regression analysis, can 
be utilized to determine the greatest similarity be- 
tween the mathematical model and the amplitude 
versus wavelength response resulting from the un- 
known, tested fluid. 

Considerable improvement in detection preci- 
sion is obtained by simultaneously utilizing at least 
several wavelengths from the entire spectral fre- 
quency range of source 17 to derive data for a 
multivariate analysis. The multivariate method al- 
lows both detection and compensation for interfer- 
ences, the detection of outlier samples, as well as 
for modelling many types of non-linearities. Since 
the calibration samples used to derive the model 
have been analyzed on a multivariate basis, the 
presence of unknown biological materials in sample 
12 does not prevent or distort the analysis. This is 
because these unknown biological materials are 
present in the pre-stored multivariate calibration 
samples used to form the model. The multivariate 
method provides for detection of outliers, i.e., those 
samples that do not exhibit the characteristic rela- 
tionship between composition and the absorbance 
spectrum of other calibration samples. Identifying 



and removing outlier samples from the calibration 
set improves the accuracy and precision of the 
analysis. Identification of the spectrum of sample 
12 as an outlier (dissimilar from the calibration set) 

5 provides methods of evaluating the validity of the 
concentration prediction calculations. Identification 
of an outlier spectrum from unknown samples, 
such as sample 12, is important since the results 
from analysis of an outlier sample are unreliable. 

to Dire consequences couid result if clinical decisions 
were to be based on unreliable analyses. It is also 
possible, utilizing known statistical methods, to 
identify the reason a sample is an outlier from the 
multivariate analysis, such as a malfunctioning ap- 

75 paratus. 

After the partial least squares technique has 
been employed by computer 39 to determine the 
characteristics of sample 12, the computer derives 
an indication of the concentration of the analyte in 

20 the unknown sample. For example, the computer 
derives the concentration of glucose in the blood- 
stream. The indication derived by computer 39 is 
coupled to conventional alphanumeric visual dis- 
play 41 . 

25 The device illustrated in Figure 1 is used to 
examine in vitro drawn samples as obtained from 
blood routed outside the body. A similar structure 
could also be used to analyze urine, saliva and 
other biological fluids for infrared active analytes 

30 thereof. 

A device operating on principles similar to 
those of the device illustrated in Figure 1, but 
which is implantable in the human body, is illus- 
trated in Figure 2. The device illustrated in Figure 2 

35 includes attenuated total reflectance crystal 51 hav- 
ing a planar face with bio-compatible crystal coat- 
ing 52 thereon. Coating 52, fabricated from any 
suitable material, such as a polyurethane foam, is 
in direct contact with the biological fluids of inter- 

40 est. Coating 52 is a porous film permitting penetra- 
tion of molecular species smaller than antibodies 
so that the molecular species of the blood or other 
biological fluid is in direct contact with the face of 
crystal 51. Alternatively, coating 52 could be thin 

45 enough to permit penetration of infrared optical 
energy from crystal 51 directly into body fluid or 
tissue. 

Crystal 51 includes inclined end faces 54 and 

55 between which coating 52 extends. End face 54 
so transmits optical energy from mid-infrared source 

56, located in housing 53 and surrounded by re- 
flector 57, having a circular cross section and 
shaped as a sheet. Reflector 57 includes mid- 
infrared source 58, adjacent end wall 54 so that 
55 infrared energy from source 56 is coupled to end 
face 54 and to the crystal interior. Infrared source 

56 is energized by DC power supply 59 in housing 
53. 
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The biological body fluid in contact with coat- 
ing 52 differentially absorbs the several 
wavelengths of wide band source 56. Thereby, the 
infrared energy at the different wavelengths of 
source 56 has variable intensity as it is propagated 
from end face 55 of crystal 51. The infrared energy 
propagating from end face 55 is generally directed 
in an upward manner, to be incident on horizontally 
extending reflective surface 61 and inclined reflec- 
tive surface 62. 

The wide band infrared energy incident on 
reflective surfaces 61 and 62 is directed to fre- 
quency dispersion device 63, thence to diode array 
detector 64, energized by power supply 59 via bus 
65. Detector array 64 includes multiple output 
leads in bus 66 so that the amplitude of the signal 
on each lead in bus 66 represents the intensity of 
the infrared energy associated with a particular 
passband associated with each detector element of 
the array. The signals in bus 66 are coupled to a 
random access memory of computer 67 which is 
constructed in the same manner as computer 39, 
Figure 1. After computer 67 determines the char- 
acteristics of analyte fluid contacting coating 52 
from the intensity versus wavelength variations in- 
cident on array 64, the computer activates an in- 
sulin pump to control the glucose concentration of 
the subject. Alternatively, the data on bus 66 could 
be transmitted to a computer external to the body 
of the subject. A wireless communication link for 
control of the pump or for the data is established 
by radio transmitter 60, connected to computer 67 
or bus 66, as required. 

Reference is now made to Figure 3 of the 
drawing wherein fiber optic element 71 , surrounded 
by a catheter (not shown) is illustrated as being 
inserted into artery 72 of arm 73 of a subject. Fiber 
optic element 71 is a single fiber optic infrared 
transmitter having an infrared reflective coating for 
internal reflection, as well known in the art. The end 
of fiber optic element 71 positioned in artery 72 
includes reflective end coating 74, while the op- 
posite end of the fiber optic element in housing 75 
has conical tip 76. 

One side of conical tip 76 is illuminated by 
wide band infrared energy either in the mid-infrared 
range or near-infrared (2.5-0.7 micrometers) range, 
as derived from a source including infrared emitter 
78 and reflector 79 having a circular sheet-like 
reflecting surface. Reflector 79 includes opening 81 
in proximity to one side of conical end 76. Optical 
energy from source 78 incident on one side of 
conical tip 76 is transmitted from the tip along the 
length of fiber optic element 71 to artery 72. The 
optical energy in fiber optic element 71 is reflected 
from end face 74 back along the length of the fiber 
optic element to conical tip 76. The portion of fiber 
optic element 71 imbedded in artery 72 does not 



include the usual protective coating such as a 
polyamide coating. Since the index of refraction of 
fiber optic element 71 is greater than that of the 
blood in artery 72 the fiber optic element functions 
5 along the entire length embedded in the artery 
basically in the same manner as attenuated total 
reflectance crystals 13 and 51. 

The infrared optical energy transmitted back 
from artery 72 via fiber optic element 71 to conical 

10 tip 76 is transmitted from the side of the conical tip 
opposite from the side adjacent opening 81 so that 
the energy transmitted from the tip is incident on 
frequency dispersion device 82, thence is coupled 
to diode detector array 83. Detector array 83 and 

75 infrared source 78 are energized by DC power 
supply 84. Detector array 83 supplies computer 85 
with signals via bus 86 to indicate the infrared 
absorption by the blood in artery 72 of the optical 
energy from source 78 for each of several fre- 

20 quency bands of the source. Computer 85 includes 
a program with a read only memory as described 
c supra with regard to computer 39, to determine the 
characteristics of the blood in artery 72. 

The optical energy from tip 76 is coupled to 

25 frequency dispersion device 82 via a reflective 
system (not shown) somewhat similar to that illus- 
trated in Figure 2. Such a reflective system is 
desirable because the emitting surface of conical 
tip 76 is off axis from frequency dispersion device 

30 82 and diode detector array 83. 

Fiber 71 has several advantages, relating to 
flexibility, ability to transmit light to remote loca- 
tions and small diameter, usually less than a mil- 
limeter, to provide intravascular, intramuscular, in- 

35 terstitial and transdermal use. With the currently 
available fiber optic elements, source 78 may be 
either a mid-infrared or a near-infrared wideband 
emitter. Currently available fiber optic elements are 
more transmissive in the 2.5-0.7 micrometer, near- 

40 infrared region than in the mid-infrared 50-2.5 mi- 
crometer region. However, the near-infrared region 
has the disadvantage of relatively low absorptivity 
by organic species relative to the absorptivity in 
the mid-infrared region. 

45 The device of Figure 3 could also be an inter- 
stitially imbedded optical fiber system or could be 
used to measure concentrations of biological fluid 
anaiytes outside of the vascular system. 

In response to the comparison made by com- 

so puter 85 of the model stored therein with the data 
derived from bus 86, the computer derives an 
indication of the characteristics of the analyte fluid 
in artery 72. The comparison and resulting con- 
centration computation are made in computer 85 

55 by the partial least squares technique. The con- 
centration indication is supplied by computer 85 to 
alphanumeric display 87. 

The model stored in computer 85 is computed 
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by a general purpose computer using the partial 
least squares method with absorption versus 
wavelength responses from plural known biological 
samples. The concentration in each of the known 
biological samples is typically determined using 
invasive enzymatic, chemical or immunological 
techniques. The model data relating absorption ver- 
sus wavelength response to concentration are read 
from the general purpose computer into the mem- 
ory of microcomputer 85 prior to installation of the 
microcomputer into the apparatus illustrated in Fig. 
3. The partial least squares computation by com- 
puter 85 also enables a determination to be made 
that the absorption versus wavelength response for 
the measured unknown sample cannot be asso- 
ciated with the concentration of the analyte in the 
blood. In particular, if the F-statistics partial least 
squares value computed from the model and un- 
known sample absorption versus wavelength re- 
sponse exceeds a predetermined value, the analyte 
concentration in the unknown sample cannot be 
determined. If this happens, display 87 is activated 
accordingly and no control over an insulin pump is 
instituted. 

In the embodiment illustrated in Figure 4, sev- 
eral infrared wavelengths in the near-infrared or 
mid-infrared range are transmitted from wide band 
source 91 in housing 92 via a fiber optic loop 
including fiber optic element 93 to ear lobe 95 of a 
human subject. The protective coating of optical 
fiber element 93 is removed to allow contact of the 
fiber with the ear as the fiber pierces through a 
portion of the ear. The exposed portion of the fiber 
surrounded by the ear thereby acts as an attenu- 
ated total reflectance element. Following penetra- 
tion of the ear, the optical fiber is again coated with 
a protective coating. The optical energy is transmit- 
ted from endface 97 to endface 101 by a continu- 
ous fiber with only region 98 extending through the 
earlobe being uncoated so there is an interaction of 
the infrared light transmitted through element 93 
with the surrounding sample medium in the ear 
lobe. 

The infrared energy emerging from end face 
101 is transmitted by horizontal and inclined planar 
reflectors 102 to frequency disperser 103 and di- 
ode detector array 104. Frequency disperser 103 
and diode detector array 104 function with com- 
puter 105, display 106 and power supply 107 in the 
same manner as described supra , with regard to 
Figs. 1-3. 

A fiber optic loop system similar to that of 
Figure 4 can be used to penetrate other portions of 
the body which permit sampling of the subdermal 
region. These types of arrangements can also be 
used with an uncoated optical fiber in a disposable 
needle device. In such a situation, the fiber is 
changed daily to prevent fibrin accumulation on the 



needle. It has been found that such a disposable 
needle penetrating the body is very well suited to 
control insulin delivery for waist-mounted insulin 
pumps. 

5 The fiber optic transmission system of Figure 4 
is more suitably used in conjunction with infrared 
sources in the near-infrared region than sources in 
the mid-infrared region. This is because the longer 
attenuation by currently available fiber optic ele- 

70 ments for mid-infrared wavelengths than for, near- 
infrared wavelengths. The devices illustrated ih.Fig- 
ures 1-4 all use attenuated total reflectance princi- 
ples. The invention is also applicable to techniques 
using transmission of infrared energy, usually in 

75 the near-infrared range. 

The present invention can be used to analyze 
biological fluids in other body parts. In the embodi- 
ment of Figure 5, used to monitor blood in the 
index finger of a human subject, housing 111 in- 

20 eludes parallel top and bottom faces 112 and 113, 
each having mating, aligned apertures 114, particu- 
larly adapted to receive a subject index finger. The 
volume in housing 111 between apertures 114 thus 
defines a cylinder, on opposite sides of which are 

25 located near-infrared source 116 and frequency 
dispersion device 119. Near-infrared source 116 is 
surrounded by sheet-like reflector 117 having a 
circular cross section. Optical energy from source 
1 16 is transmitted directly toward frequency disper- 

30 sion device 119 and is reflected by reflector 117 so 
that rays of the near-infrared energy are transmit- 
ted through an index finger inserted into housing 
111. The rays transmitted through the finger are 
incident on frequency dispersion device 119 and 

35 on parallel horizontally disposed reflecting surfaces 

120 which directs them to device 119. In close 
proximity to frequency dispersion device 119 is 
diode detector array 118 which supplies computer 

121 with signals to enable display 122 to be en- 
40 ergized. Power supply 123 energizes the diodes of 

array 118 and lamp 116. 

Hence, in the embodiment of Figure 5, optical 
energy in the near-infrared spectrum is transmitted 
directly through the index finger inserted through 

45 apertures 114 into housing 111. Blood in the index 
finger differentially absorbs different wavelengths 
from source 116. The differential absorption is de- 
tected by monitoring the amplitude at the different 
wavelengths coupled by frequency dispersion de- 

50 vice 1 19 to diode detector array 118. 

Principles of the direct near-infrared transmis- 
sion system illustrated in Figure 5 can be used to 
measure concentrations of species in other body 
parts. For example, housing 111 can be modified 

55 so that the abdomen, ear lobe, wrist or other ap- 
pendages are located between the near-infrared 
source and reflector combination and the frequency 
dispersion device and diode detector array. 
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The device illustrated in Figure 6, used to 
monitor biological fluids non-invasively using near- 
infrared spectroscopy and diffuse reflectance sam- 
pling, includes housing 131 in which are located 
near-infrared source 132 and reflector assemblies 5 
133 and 134. Reflector assembly 133 is configured 
as a sheet having a circular cross section that 
approximately surrounds source 132 and has an 
opening at one end face of housing 131. Reflector 
assembly 134, located immediately below reflector io 
assembly 133, includes an opening in the same 
end face of housing 131. The end face of housing 
131 where the apertures of reflector assemblies 
133 and 134 are located is adapted to be placed 
directly against a body part, such as the head or is 
the abdomen, for monitoring biological fluids in the 
brain or liver. The optical energy from source 132 
and reflector assembly 133 penetrates to the inte- 
rior of the body part of interest and is diffusely 
reflected from the blood or other biological fluid in 20 
the organ being monitored. The optical energy 
from source 132 and reflector 133 is differentially 
absorbed by the organ being analyzed within the 
body, so that diffusely reflected optical energy 
propagated from the organ through the body is 25 
incident on reflector assembly 134. The diffusely 
reflected infrared energy incident on reflector as- 
sembly 134 is coupled to frequency dispersion 
device 135, thence to diode detector array 136 
which drives computer 137 to energize display 138, 30 
as described supra . Power supply 139 energizes 
infrared source 132 and diode detector array 136. 

Reference is now made to Figure 7 of the 
drawing, a schematic diagram of an implantable 
apparatus for monitoring biological fluids using 35 
near-infrared spectroscopy, wherein near-infrared 
energy is transmitted through an internal body part, 
such as artery 141. The apparatus illustrated in 
Figure 7 is located in housing 142, including par- 
allel top and bottom faces 143 and 144 having 40 
mating circular apertures through which artery 141 
is adapted to fit. Housing 142 includes hinged end 
portion 245. connected to the remaining, main por- 
tion 148 of the housing by hinges 246. Housing 
portions 245 and 148 include a latch arrangement 45 
(not shown) so that artery 141 can be surrounded 
by the housing and fit into the mating circular 
apertures. 

End portion 245 of housing 142 includes near- 
infrared source 146, surrounded by sheet-like re- 50 
flective element 147, having a circular cross sec- 
tion. Main portion 148 of housing 142 includes 
parallel planar, horizontally disposed reflective ele- 
ments 149 and 150 which direct optical energy 
transmitted from source 146 and reflector 147 55 
through artery 141 onto frequency dispersion de- 
vice 152. Frequency dispersion device 152 is lo- 
cated in the main part of housing 142 and is 



optically coupled to diode detector array 1 53 which 
supplies signals to computer 154 that communi- 
cates with an insulin pump via a wireless link 
including radio transmitter 250. Alternatively, the 
signals from diode detector array 153 are coupled 
to an external computer via the wireless link. Near- 
infrared source 146 and the diodes of detector 
array 153 are energized by DC power supply 155. 

Devices similar to that illustrated in Figure 7 
could be used for transmission through other tis- 
sues, such as the liver, for examining bile, or the 
bladder for urine analysis. A further variation of the 
implanted device of Figure 7 involves the use of an 
external near-infrared source, to conserve energy 
of the power supply. In such an instance, optical 
energy from the near-infrared source is transmitted 
through a tissue sample to an implanted detection 
apparatus. Infrared energy transmitted through the 
tissue is incident on a frequency dispersion device 
and a diode detector array. 

Three preferred embodiments of the frequency 
dispersion device and the diode detector array of 
any of Figures 1-7 are illustrated in Figures 8-10. In 
each of the embodiments of Figures 8-10, infrared 
energy is incident on a passive elongated diode of 
the detector array. Each diode has associated 
therewith a certain wavelength range, either in the 
mid-infrared or near-infrared spectrum. The inten- 
sity of the infrared radiation incident on each semi- 
conductor diode changes the charged carrier den- 
sity of the diode to control the impedance thereof 
between the diode anode and cathode electrodes. 
Each diode is energized by a DC voltage from the 
power supplies of Figures 1-7, so that the voltage 
output of each diode varies as a function of the 
intensity of the incident infrared energy in the band 
of interest for the particular diode. 

To these ends, diode detector array 200. Fig- 
ure 8, includes several elongated, mutually insu- 
lated semiconductor diodes 201-211, having mutual 
planar faces on which is superimposed filter frame 
213 including an equal number of several individual 
elongated infrared filter elements 221-231. Each of 
filter elements 221-231 is superimposed on and 
abuts against a corresponding one of diodes 201- 
211. Each of elongated filters 221-231 has a pass- 
band for a different wavelength interval in the mid- 
infrared or near-infrared spectrum. Thereby, a volt- 
age is derived across the electrodes of each of 
diodes 201-211, which voltage is indicative of the 
intensity of the radiation in the infrared region 
passed by the filter abutting against the particular 
diode. The signals across the electrodes of diodes 
201-212 are supplied to a bus, such as bus 38, 
Figure 1. which couples the signals to designated 
locations in the computer random access memory. 

A similar result is obtained with the structure of 
Figure 9, wherein diode detector array 300 is iden- 
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tical to array 200, Figure 8, and thereby includes 
several parallel, elongated semi-conductor diode 
elements 301-311. Superimposed on the elongated 
planar faces of the diodes of array 300 is planar 
sheet 313, fabricated of a variable frequency re- 
sponse material. The material of sheet 313 is ar- 
ranged so that the material in the sheet on one 
side (e.g., the left side thereof, as illustrated in 
Figure 9), is transparent only to the shortest 
wavelength of interest in the near-infrared spectrum 
emitted by the near-infrared source or only the 
shortest wavelength of interest in the mid-infrared 
spectrum emitted by the mid-infrared source. At 
the right side of sheet 313, the material of the 
sheet is transparent only to the longest 
wavelengths of interest in the near-infrared or mid- 
infrared spectra of the source. 

The intermediate portions of sheet 313 are 
fabricated of material that is transparent only to 
intermediate portions of the near-infrared or mid- 
infrared spectra. The passbands of the material in 
sheet 313 progressively increase in wavelength 
from the left to the right edge of the sheet. Sheet 
313 abuts against elements 301-311 of array 300 
so that diode 301 is responsive basically to the 
same range of wavelengths as diode 201 in array 
200 of Figure 8. Similarly, elements 302-311 of 
array 300 derive responses for the same 
wavelengths as corresponding elements 202-21 1 of 
array 200. 

Similar results are attained with the detector 
array of Figure 10, wherein array 400 includes 
several elongated semiconductor detectors 401- 
411. Each of the diodes 401-412 is differentially 
doped so that it is responsive to a different 
wavelength in the near-infrared or mid-infrared 
spectrum. Hence, diode 401 is fabricated or doped 
so that it is responsive only to infrared energy 
having the shortest wavelength of interest in the 
near-infrared spectrum or in the mid-infrared spec- 
trum. Diode 411 is doped so that it is responsive 
only to the longest wavelengths of interest in the 
near-infrared or mid-infrared spectra. Intermediate 
diodes 402-410 are fabricated or doped so that 
they are responsive to successively longer 
wavelengths in the near-infrared or mid-infrared 
spectra. 

The multivariate approach of the present inven- 
tion permits the concentration of molecular sub- 
stances in fluids to be predicted without the use of 
enzymatic, chemical or immunological methods. 
The partial least squares technique produces re- 
sults as precise as reactive determinations by the 
enzymatic, chemical or immunological methods. By 
using the multivariate mid-infrared and near-infrar- 
ed techniques, it is possible with the invention to 
determine the concentrations of substances without 
penetrating the body, except for the optical radi- 



ation. Such a non-invasive determination has a well 
defined need in the control of diabetes mellitus. 
The invention has use for determining the levels of 
alcohol, ketones, fatty acids, cholesterol, 

5 lipoproteins, triglycerides, white blood cells, al- 
bumin, blood urea nitrogen, creatinine, concurrent 
medications, such as drugs, inorganic molecules 
such as phosphates, and detection of other infrared 
active compounds. The disclosed configurations 

to can be continuously or intermittently operated to 
transmit concentration or dosage data to an im- 
planted or external pump for drug delivery, such as 
for insulin. 

While there have been described and Mus- 
ts trated several specific embodiments of the inven- 
tion, it will be clear that variations in the details of 
the embodiments specifically illustrated and de- 
scribed may be made without departing from the 
true spirit and scope of the invention as defined in 
20 the appended claims. 



Claims 

25 1. A method of determining one or more un- 
known values of known characteristics, such as the 
concentration of an analyte, of a biological fluid in a 
mammal, said method comprising: 

(a) irradiating said biological fluid having said 
30 unknown values of said known characteristics with 

infrared energy having at least several wavelengths 
so that there is differential absorption of at least 
some of said wavelengths by said biological fluid 
as a function of said wavelengths and said char- 
as acteristics, said differential absorption causing in- 
tensity variations of said wavelengths incident from 
said biological fluid as a function of said 
wavelengths and said unknown values of said 
known characteristics, or as a function of said 
40 wavelengths and said characteristics, or as a func- 
tion of said characteristics having unknown values; 

(b) measuring said intensity variations from 
said biological fluids; and 

(c) calculating said unknown values of said 
46 known characteristics in said biological fluid from 

said measured intensity variations utilizing an al- 
gorithm and a mathematical calibration model, said 
algorithm including all independent sources of in- 
tensity variations v. wavelength information ob- 

50 tained from irradiating a set of samples in which 
said values of said known characteristics are 
known, said algorithm also including more 
wavelengths than samples in said set of samples, 
said model being constructed from said set of 

55 samples and being a function of said known values 
and characteristics and said intensity variations v. 
wavelength information obtained from irradiating 
said set of samples. 
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2. The method of claim 1 wherein the fluid is 
irradiated in vivo and non-invasively. 

3. The method of claim 1 wherein the fluid is 
irradiated in viva and invasively. 

4. The method of claim 1 wherein the fluid is 5 
irradiated in vitro. 

5. The method in any of claims 1-4, wherein 
said infrared energy is near infrared energy. 

6. The method of any of claims 1-5, wherein 

said algorithm is selected from the group including, 10 
partial least squares and principle component re- 
gression. 

7. The method of either of claims 2 or 3, 
wherein said biological fluid is irradiated by direct- 
ing said infrared energy so that it is incident on a 75 
portion of said mammal which includes said bio- 
logical fluid, so that said biological fluid partially 
absorbs said incident infrared energy. 

8. The method as set forth in claim 7. wherein 

said infrared energy incident on said portion of said 20 
mammal passes through said portion. 

9. The method as set forth in claim 8, wherein 
said infrared energy incident on said portion of said 
mammal is partially absorbed by said portion and 
partially diffusely reflected from said portion. 25 

10. The method as set forth in any of claims 1- 
9, further including the step of detecting outlier 
samples in said set of samples. 

11. The method as set forth in claim 10, 
wherein said detection of outlier samples in said 30 
set includes the step of comparing said intensity v. 
wavelength responses of each sample n said set to 

said model, and determining a measure of the 
magnitude of the difference between said intensity 
v. wavelength responses of each of said samples 35 
and said model in response to the comparing step. 

12. The method as set forth in claim 11, 
wherein said detection of outlier samples in said 
set includes comparing said known values of said 
characteristics of each sample in said set to the 40 
calculated value of said characteristic, said cal- 
culated value being based on said model estimate 

of said value of said characteristic for each of said 
samples, and determining a measure of the mag- 
nitude of the difference between said known value 45 
of said characteristic and said calculated value of 
said characteristic. 

13. The method as set forth in either of claims 
11 or 12, further comprising indicating the probabil- 
ity of said magnitude being caused by random 50 
chance by using a statistical test, and classifying 

as outliers said samples in said set having an 
excessively low determined probability. 

14. The method as set forth in claims 12 or 13 
wherein the statistical test is the F ratio test. 55 

15. Apparatus for determining at least one un- 
known value of a known characteristic, such as the 
concentration of an analyte, in a biological fluid. 



said apparatus comprising: 

(A) a source of infrared energy having at 
least several wavelengths 

(B) means for coupling said at least several 
wavelengths of said infrared energy to said biologi- 
cal fluid to enabie said biological fluid to differen- 
tially absorb at least some of said wavelengths, 
said differential absorption causing intensity vari- 
ations of said infrared energy incident from said 
biological fluid as a function of said at least several 
wavelengths of said energy and said unknown val- 
ue of said known characteristics; 

(C) means for measuring said intensity vari- 
ations; and 

(D) computer means including: 

i. a stored model constructed from a set of 
samples in which the values of said known char- 
acteristic are known, said model being a function of 
said known values from said set of samples and 
intensity v. wavelength information obtained from 
said set of samples, and 

ii. an algorithm for utilizing (a) all indepen- 
dent sources of said intensity variations v. said 
wavelengths information from both said set of sam- 
ples and said biological fluid and (b) more 
wavelengths than samples, said algorithm utilizing 
said model for calculating said unknown value of 
said known characteristic of said biological fluid 
from said measured intensity variations from said 
biological fluid. 

16. The apparatus as set forth in claim 15, 
further including means for determining whether 
said intensity variations v. wavelength of said 
known characteristic in said biological fluid is an 
outlier. 

17. The apparatus as set forth in either claim 
15 or 16, wherein said means for coupling includes 
an attenuated total reflectance (ATR) device. 

18. The apparatus as set forth in claim 17, 
further including a flow cell, said ATR device being 
positioned in said cell for in vitro sampling. 

19. The apparatus as set forth in claim 17, 
wherein said ATR crystal has a biocompatible coat- 
ing on at least a portion thereof whereby said 
portion of said ATR device can be placed in con- 
tact with said biological fluid in vivo. 

20. The apparatus as set forth in either claim 
15 or 16, wherein said means for coupling includes 
a fiber optic device. 

21. The apparatus as set forth in claim 20, 
wherein a portion, of said fiber optic device forms 
an ATR device. 

22. The apparatus as set forth in claim 20. 
wherein said fiber optic device includes a portion 
for transmitting said infrared energy to said biologi- 
cal fluid and a separate portion for transmitting said 
infrared energy from said biological fluid to said 
apparatus. 
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23. The apparatus as set forth in either of 
claims 15 or 16, wherein said apparatus includes a 
first portion which can be positioned on one side of 
an in vivo source of biological fluid and a second 
portion which can be positioned on another side of 5 
said in vivo source of biological fluid. 

24. The apparatus as set forth in either of 
claims 15 or 16, wherein said means for coupling 
includes means for transmitting said infrared en- 
ergy to an in vivo source of said biological fluid to 
and means for measuring diffuse infrared reflection 
from said in vivo source. 

25. The apparatus as set forth in either of 
claims 15 or 16, further including means for trans- 
mitting signals from said computer means to a 75 
means for changing said value of said known char- 
acteristic. 

26. The apparatus as set forth in claim 25, 
wherein said means for changing is an insulin 
pump. 20 

27. The apparatus of claim 15 or 16, wherein 
said means for measuring said intensity variations 
includes an array of infrared sensors, means for 
frequency dispersing said intensity variations of 

said at least several wavelengths onto said sen- 25 
sors, different ones of said sensors being provided 
for different ones of said at least several 
wavelengths, so that said different ones of said 
sensors derive different electric signals. 
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